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Abstract 

We calculate parity nonconserving observables in the processes where a neutron is captured on 
a proton at the threshold energy radiating a photon. Various potential models such as Paris, 
Bonn and Argonne vl8 are used for the strong interactions, and the meson-exchange description is 
employed for the weak interactions between hadrons. The photon polarization P 7 in the unpolarized 
neutron capture process and photon asymmetry Ay in the polarized neutron capture process are 
obtained in terms of the weak meson- nucleon coupling constants. turns out to be basically 
insensitive to the employed strong interaction models and thus can be uniquely determined in terms 
of the weak coupling constants, but P 7 depends significantly on the strong interaction models. 
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I. INTRODUCTION 



Recent experiments to explore the weak interactions between hadrons through parity 
nonconserving (PNC) observables in nuclear systems or reactions [3|, |4j have triggered 
revived interests in this field. These PNC observables can be related to the meson-nucleon 
weak coupling constants which are introduced in the meson-exchange potential description 
of the hadronic weak interaction 5]. However, due to various uncertainties (see Ref. (s| 
for details), the weak coupling constants were fixed only within certain ranges [5j. Thus, it 
has been hoped that the PNC observables from various experiments can reduce the range of 
these coupling constants and eventually determine the values. The situation, however, has 
not been much improved even by the recent measurements. For example, the value of the 
7i — N weak coupling constant, h\, from the anapole moment of 133 Cs |l| is inconsistent with 
a previous value obtained from the forbidden 7-decay of 18 F (|: h\ determined by the former 
is larger than that from the latter by a factor of 7. New experiments, already completed (^J, 
being done (4j or expected to be performed, concern two-nucleon systems in which many 
body effects are absent. Thus they are expected to give more stringent constraints on the 
weak coupling constants. For the current status of research on the weak coupling constants, 
see 

In this work, we calculate the photon asymmetry A 1 in the radiative capture of a polarized 
neutron on a proton, n+p — > cZ+7, and the circular polarization of photons P 1 in n+p — > d+7 
at the threshold. The latest experimental value of A 1 is —(1.5 ± 4.8) x 1CT 8 j8j, but the 
experiment being done at LANSCE aims at the accuracy of 10~ 9 0]. Theoretical calculations 
of A 1 using strong models made in the 1960's and the 1970's such as Hamada- Johnston, Reid- 
soft-core and Tourreil-Sprung show results similar to each other; v4 7 ~ — O.ll/i* jsj. Ay is 
predominantly determined by h\ and depends very little on other coupling constants (as 
will be shown in Table H}. In this work we present A 1 calculated with potentials such as 
Paris Q, Bonn Bonn-A and Bonn-B |l2|, and Argonne t>18 (Ai>18) |l^J. We compare 
our results with previous ones Q and investigate the model dependence of A 1 . 

Contrary to Ay, P~ at the threshold is known to be sensitive to the heavy meson (p and 

on 

u>) components of the weak potentials |9j,llJ|- The most recent experimental value of P 7 is 
(1.8 ± 1.8) x 10 -7 and theoretical calculations made in the 1970's agree with this value 
within the experimental errors. However, since P 7 is sensitive to the short range properties 
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of the strong interactions as well as of the weak interactions, its model dependence is more 
noticeable than A~ Since the inverse process, 7 + d — > n+p, whose PNC asymmetry 

at the threshold is equal to P 7 , becomes experimentally feasible nowadays, we expect that 
P 7 can be measured more precisely and can provide more constraints on the weak dynamics 
of hadrons. We thus investigate the model dependence of P 7 with the same potentials that 
we use in calculating A T 

In Sect. 2, we present the Desplanques-Donaghue-Holstein (DDH) potential [5] and the 
parity-admixed wave functions in the initial and the final states. In Sect. 3, the electromag- 
netic operators are presented, matrix elements are derived, and the results for and P 7 
are shown. Discussions on the results follow in Sect. 4. 

II. PARITY ADMIXED WAVE FUNCTION 

The Schrodinger equation for a two-nucleon system can be written as 

1 (id 2 i(i + iy 



Hmr) 



r - n I + V c {r) + V T (r) S 12 (r) + V pnc (r 



*(r) 



rriN \r dr 2 r 2 
E*(r), (1) 



where Vr represents the tensor potential and Vc includes central, spin-orbit, spin-spin and 
quadratic spin-orbit interactions in the strong potential. In the Paris and Bonn potentials, it 
is essential to include the momentum dependent term in the central potential to obtain the 
correct phase shifts even at low energies. In Ref. a transformation useful for treating 
the momentum dependent term is suggested. In this work, however, we have dealt with 
the momentum dependent terms without using such a transformation and have solved the 
Schrodinger equation as it is. We have confirmed that the solutions thus obtained reproduce 

nnn 

the results of each potential model [I1J, llJj, Il2j fairly well with differences less than 1%. The 



small differences can be attributed to the use of slightly different values of physical quantities 
in the calculations. V pnc is the PNC potential, and we use the one given by DDH [5] 

V pnc (r) = V; nc (r) + V; nc (r) + V p " nc (r), 

V p Ur) = i|=^ (n x r 2 f (<t x + <r 2 ) • [p, f^)] , (2) 



VP (r) = 9pNN 
pnA ' ~" m N 



(h° p n • r 2 + \h\ (7f + r|) + i (3rM - r x ■ r 2 ) 
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((<ri - <r 2 ) • {p, f P (r)} + i(l + *„) (<T! x 0-2) • [p, f p (r)} ) 
+ ^ (ri x t 2 ) z {a, + <r 2 ) • [p, /,(r)] - ^ (r* - r 2 *)(«n + <r 2 ) • {p, f p (r)} 



,(3) 



V" (r) 

' pnc \ J 



9wNN_ 

m N 



{hi + \hl (7f + r|)) x 



((m - <r 2 ) • {p, / w (r)} + i{\ + X .) (*i x <r 2 ) • [p, / w (r)] ) 
+^W-^)(*i + ^)-{p,/ w (r)} 



(4) 



where the strong coupling constants are g n NN = 13.45, q p nn = 2.79, g U NN = 8.37 and the 
anomalous magnetic moments are Xp — 3-71 and Xu — —0.12. The Yukawa functions /m(?") 
are defined as 

e ~m M r 



f M (r) = 



(M = tt, p, u). 



Anr 

The quantities h^ 1 represent the weak meson-nucleon coupling constants where AI denotes 
the isospin transfer. 

At the threshold energy, the initial scattering state, n + p, is dominated by the lowest 
angular momentum state, i.e., the 1 S channel, and higher angular momentum states are 
suppressed. Thus in this work we just include the next low-lying state, the 3 Si — 3 D± partial 
waves, where the 3 D\ state is induced by the tensor interaction in the initial scattering 
state. Then the parity-even state of the initial wave function consists of the 1 S , 3 Si and 
3 Di states. 

Since V pnc is a parity-odd operator, it creates opposite parity components in the wave 
function. For example, when V pnc is operated on the 1 S state, the isoscalar and isotensor 
terms of V pnc generate a 3 Po admixture, where the tilde denotes the parity-admixed com- 
ponents generated from the DDH potential. Similarly, 3 Pi and 1 P 1 admixtures arise from 
applying the isovector and isoscalar components of V pnc to the 3 Si — 3 Di state, respectively. 
The total wave function of the initial state with its parity admixture at the threshold can 
be written as 



* 4 (r) = *f(r) + *r c (r 
1 



mr) = 



4:71V 
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V8 



wt{r) Xis z C 



00 



Aixr 



- v?\r){<r x + <r 2 ) X is z Go + \ vf\r){cr x - <r 2 ) X oo C 



10 



(5) 
r, (6) 
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where Xss z and (tt z represent the spin and isospin part, respectively. u s is the radial part 
of the wave function for the 1 5'o channel, Ut for 3 Si and Wt for 3 D\. The final state wave 
function can be written in a similar way as 



/I 
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(7) 
r, (8) 



where «d(r) (it^r)) is the radial wave function for the 3 Si ( 3 -Di) deuteron state, and v/ and 
Vj pl denote the parity nonconserving admixture due to the 1 P\ and 3 P\ states, respectively. 
By inserting the initial and the final wave functions into the Schrodinger equation ((TJ with 
the strong and weak PNC potentials, one can obtain the radial wave equation for each 
channel (see Appendix for details). 



III. MATRIX ELEMENTS, P 7 AND A 7 

At the threshold energy, it is well known that the neutron capture cross section is domi- 
nated by the isovector Ml transition. We can evaluate the parity conserving Ml transition 
amplitude by using the one-body spin current operator 



' M 



X>f<TiXfc 7 , (9) 



4m A? 

where /j,y = 4.71, and fc 7 is the photon momentum. Amplitudes between the states with 
opposite parities would become non-zero through the El transition. While the impulse 
approximation is used in evaluating the Ml amplitude, the contribution from the exchange 
currents can be well accounted for by the Siegert's theorem. The El current operator with 
Siegert's theorem reads 

J§ = -ij W - )r, (10) 
where u is the photon energy (2.2246 MeV at threshold). The transition amplitudes (M-{) 
can be classified in terms of the electromagnetic type Ai(=M or E), the initial (i) and the 
final (/) states. The leading parity-conserving isovector Ml transition occurs between the 
initial 1 S'o and final 3 Si states, and we denote its amplitude by Mf^. The non-zero PNC 
El amplitudes are represented similarly as E^ for 3 P — > 3 S\ + 3 D 1 , E^ Q for 1 Sq — ► l P%, 



El°\ for 3 Pi -> 3 S'i + 3 £>i, and P^ 1 for 3 5i + 3 P>i -> 3 Pi, where the tildes are to distinguish 
the PNC amplitudes from the normal parity conserving ones. With the wave functions of 
Eqs. (JSJ)-©, we obtain the matrix elements 

Mt'o = ^ - [dru d (r)u s (r), (11) 



N 



E 3 4l = drr(u d (r) - V2w d (r)) «*°(r), (12) 



#1$ = *r^ pl (r)u s (r), (13) 

In terms of these electromagnetic amplitudes, the two PNC observables are written as 

a i = -2 *;If 3 " ^ 4 + 4, us) 

Jw lsO 

p , _ ^3 P + ^UO _ pi p f , . 

m lsO 

where A\ = -2Ef p \ / Mf^ and P\ = -2E^l/M?^ have the PNC component in the initial 
state and A? = -2E^ s l/Mf s J and P* = -2E%& / 'Mf s J have the PNC component in the final 
state. Numerical results are given in Tabled We express the results for A^ and P 7 in terms 
of the weak coupling constants hj^ to show explicitly the dependence of A 7 and P 7 on each 
meson. "Best values" refer to A^ and P 7 values evaluated with the so-called best values of 
the weak meson-coupling constants suggested by Ref. p. They are h? = —11.4, h!^ = —1.9 
hi = -9.5, hi = 4.6, h\ = -0.2 and hi = -1.1, in units of 1(T 7 . 



IV. RESULTS AND DISCUSSIONS 

Asymmetry (A 7 ) 

As shown in Table HI the Bonn and Avl8 models predict the same v4 7 value, while the 
best value from the Paris potential is larger in magnitude than those from Bonn and Aw 18 
by a factor of 1.27. This factor can be understood by examining the wave functions that 
contribute to A^. u s and u d are plotted in Fig. HI and v^ pl in Fig. El As can be seen in Fig.[T] 
and El u d and v^ pl calculated with different potentials are very similar to each other for all 
potentials, but u s calculated with the Paris potential is substantially different from u s from 
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Model 


I 


Best values (xlO -8 ) 




4 


4 

/ 


Paris 


-0.148/4 - O.OOl/i* + 0.003/ii, 


-6.85 


-3.34 


-3.51 


Bonn 


-0.117hl - O.OOl/i* + 0.003/ii, 


-5.42 


-2.66 


-2.76 


Bonn-B 


-0.117/4 - O.OOl/i* + 0.002/ii, 


-5.41 


-2.65 


-2.76 


Avl8 


-0.117hl - O.OOl/i* + 0.002/ii, 


-5.41 


-2.63 


-2.78 


Model 


p 7 


Best values (xl0~ 8 ) 


P, 


pi 

7 


Pf 

7 


Paris 


-0.0106/1° + 0.0074/1° - 0.0191/i2 

fj U/ fj 


2.88 


-1.24 


4.12 


Bonn 


-0.0890/i° + 0.0088/i° - 0.0214/1* 


12.01 


-1.40 


13.4 


Bonn-B 


-0.0286/i° + 0.0012/i° - 0.0208/^ 


5.21 


-1.35 


6.56 


Avl8 


-0.0088/i° + 0.0034/i° - 0.0175/^ 


2.64 


-1.11 


3.75 



TABLE I: Results for the observables and P 7 for various phenomenological models in terms 
of the weak coupling constants h^J . Best values mean Ay and P 7 values obtained with the best 
values of suggested by DDH. They are in units of 10~ 8 . The definitions of P^ and A 1 <1 are 
given in Eqs. (|17j) and (j!6j) . respectively. 




r (fm) r (fm) r (fm) 

FIG. 1: The wave functions u s (r), n^(r) and Wd(r) calculated with different potentials are plotted. 
The results for u s (r) from Bonn, Bonn-B and At>18 are indistinguishable and correspond to the 
upper curve in the figure. ut(r) and wt(r) are not shown here, because they do not depend very 
much on the models. 
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r (fm) r (fm) r (fm) 



FIG. 2: Wave functions for the parity- admixed states for different potentials. The wave functions 
are given in units of h~. Note the difference in the scale. 

other potentials. {vl pl , though not shown here, is more or less the same for all potential 
models.) The reason for this difference can be traced back to the fact that the Paris potential 
was fitted to proton-proton data and therefore yields a scattering length of a = —17.6 fm 
for the 1 S'o channel (in the absence of the Coulomb interaction) while Avl8 and Bonn are 
neutron-proton models and yield a = —23.7 fm. The Ml isovector amplitude Mf^ in the 
denominator of Eqs. ()16j) and (JT7jl is 0.184 for Paris potential and 0.233 for Av 18 and Bonn 
potentials. Since the wave functions that contribute to the El amplitude in Eqs. ()14j) and 
f!15|) (numerator of Ay) are very similar, the difference in the values of Ay comes mostly 
from the value of the Ml amplitude in the denominator of Ay. Indeed, the ratio of the Ml 
amplitudes for Awl8 to Paris, 1.27, is equal to the ratio of the best value of Ay. Thus, if 
one could readjust the Paris potential to produce the accepted n — p scattering length in 
the 1 So channel, the four models would give us essentially model-independent Ay values. A 
recent work Il8f in which various contributions from the exchange currents are taken into 

n 

account confirms this model-independent nature of Ay. The magnitude of Ay in Ref. [18|| 
(Ay = —4.98 x 10 -8 ) with pion-exchange currents is smaller than ours (Ay = —5.41 x 10~ 8 ) 
by about 9 %, and is in agreement with the result of Ref. (Ay = —4.94 x 10 -8 ) where 
one-body and leading pion-exchange currents are considered. On the other hand, if one 
employs the h\ value from the 18 F |6j and 133 Cs [lj experiments, Ay becomes —1.52 x 10~ 8 
and —11.1 x 10~ 8 , respectively. Since the contribution from the pion to Ay is more than 
99 % of the total value (see Table HJ), a more accurate measurement of Ay can provide a 
stringent determination of h\. 
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FIG. 3: Integrands that enter into the El amplitude of P 7 . The left panel shows I\(r) = r (ud(r) — 
\f2wd{r)) v 3p0 (r) of Eq. (|12|) . and the right one l2{r) = rv^ pl (r) u s (r) of Eq. 

In passing, we remark that A 1 from a previous work j^J using Paris potential differs from 
our A 7 in sign though the magnitudes agree. It appears that the definitions of differ in 
sign. 

Polarization (P 7 ) 

While A 7 is dominated by the long range part of the interactions and is practically model- 
independent, P 7 depends strongly on the heavy meson exchanges and on the potential model. 
P 7 's calculated with the best values of the weak coupling constants [5] and Paris and Av 18 
potentials are similar to each other, but P 7 's evaluated with Bonn and Bonn-B are larger 
than that with Avl8 by a factor of 5 and 2, respectively. (Bonn-A and Bonn produce similar 
results and thus Bonn-A is not included in the discussion.) P 7 's expressed in terms of fiff in 
Table U show that P 7 from Bonn is more sensitive to h° p than P 7 from other potentials, while 
the terms depending on h° u and h 2 are only moderately model dependent. The contributions 
from the initial ( 3 Po) and the final ( 1 Pi) states listed in Table |U show that the initial state 
contribution (P*) is rather model- independent, but the contribution from the final state 
{Pi) is highly dependent on the potentials. 

The numerical factors in front of the weak coupling constants in Table |U are determined 
by the strong potentials in each channel through the wave functions of the 1 Sq and 3 Si — 3 Di 
channels that enter into the source terms in the Schrodinger equation of the parity-admixed 
states (see Appendix). The wave functions for the 3 Po and the 1 P\ states are shown in 
Fig. El and the corresponding integrands of the El amplitudes, Eqs. (|T2"|) and (fT3|). are 
shown in Fig. El The wave function (S 3p0 (r)) and the integrand (I\(r)) for the 3 Po state 
exhibit a sizeable model dependence. However, for the x Pi channel there are even more 
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FIG. 4: Comparison of the central potentials in the P\ channel of the Bonn, Bonn-B, and Av 18 
potentials. 
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FIG. 5: Source terms for the 1 Pi and 3 Pi states. Compared to the terms of 3 Pi, those of l P\ are 
more model dependent. Note the difference in the scale. 

drastic variations in both wave function (v^ pl (r)) and the integrand (I 2 (r)) depending on 
different potentials. Such a strong model dependency can be understood from the behavior 
of the strong potential in the 1 P\ channel and the source term that contributes to v 1 / 1 in the 
Schrodinger equation (see Eq. (|A.6)0 . Fig.|3]shows the strong potentials in the l P\ channel. 
The Bonn potential for l P\ channel becomes attractive in the short range region while Aw 18 
is repulsive in the whole region. The attraction at short ranges increases the probability for 
a nucleon to be present in the region, and this can partly explain the shape of v 1 / 1 in Fig. |21 
A recent work by R. Schiavilla et al. shows a similar behavior of v 1 / 1 . 

In Fig.0 we compare the source terms of the 3 Pi (the right hand side of Eq. (jA.5|) ) and 1 P± 
(the right hand side of Eq. (jA.6|) ) states. The sources for the 3 Pi state from different models 
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show a moderate model dependence for r < 0.5 fm, but they have significant magnitudes and 
are indistinguishable in the intermediate and long range regions, which explains the model 
independent results of in Table HI On the contrary, most of the contribution to the source 
terms of v d pl comes from the intermediate and short range region, and they depend strongly 
on the model. A larger source combined with attraction in the short-range region, as it is 
the case for the 1 P\ channel of the Bonn models, yields an enhanced contribution to P T 

Concluding, we have calculated parity non-conserving observables P 7 for the reaction 
n + p — > d + 7 and for the reaction n + p — » cf + 7 at threshold. We have employed 
the Paris, Bonn and Av 18 potentials for the strong interaction and the DDH potential for 
the weak interaction. A 1 turns out to be independent of the strong interaction models, 
while P 7 is sensitive to the dynamics at short ranges. Since v4 7 is rather strong-interaction 
independent, one can reduce the uncertainty in the value of h\ by measuring v4 7 accurately. 
Regarding P 7 , there are relatively large uncertainties, which stem from ambiguities in both 
strong and weak interactions at short ranges. However, since the major uncertainty comes 
from the X P\ channel and the value of h° p , an accurate experimental measurement of P 7 can 
shed some light on the weak coupling constants. 
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APPENDIX 

The radial equations for the 1 Sq continuum, and the 3 Si — 3 Di continuum and bound 
states read 



u" s {r) + m N (E- V c (r)) u s (r) 
u t(d)( r ) + m N (E- V c (r)) u t{d) {r) 
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w 't(d)( r ) - -^ w t{d)(r) -m N (E- V c {r) + 2V T {r))w t{d) {r) 
The equations for the parity-admixed states are 



vf P °"(r) - - 9 vr + m N (E- V c {r) + W T (r)) v\ 



~3p0 



~3p0 



0, (A.l) 

y/8m N V T (r)w t{d) (r),(A.2) 
V8m N V T {r)u t(d) {r). (A.3) 



11 



-2 ( Xp + 2) u s {r)^ (f° p (t) - \/|^ 2 (r) j + (x. + 2) «,(r) ^(r) 
+2r (w-^(r) + ^(r))A(!^) 



2 

73 



U t (d)( 



+2v / 2(F p 1 ( 



2^2 



(^(0-^(0) (u t{d) (r) - V2w t{d) (', 



vr\r)-^vr + rn N {E-V c {r))vr = 



2 

7! 



(u d (r) - V2w d (r)) (3 Xp F° p 
-2 (3F p °(r) - F°(r)) |- (u d (r) - V2w d (r)) 



where F^(r) = g MNN h^f f M (r). 



(A.4) 



(A.5) 



(A.6) 
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